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We apply the approach of the Google matrix, used in computer science and World Wide Web, to 
description of properties of neuronal networks. The Google matrix G is constructed on the basis of 
neuronal network of a brain model discussed in PNAS 105, 3593 (2008). We show that the spectrum 
of eigenvalues of G has a gapless structure with long living relaxation modes. The PageRank of the 
network becomes delocalized for certain values of the Google damping factor a. The properties of 
other eigenstates are also analyzed. We discuss further parallels and similarities between the World 
Wide Web and neuronal networks. 
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I INTRODUCTION 

More than 50 years ago John von Neumann traced first 
parallels between architecture of the computer and the 
brain [1] . Since that time computers became an unavoid- 
able element of the modern society forming a computer 
network connected by the World Wide Web (WWW). 
The WWW demonstrates a continuous growth approach- 
ing to 10 11 web pages spread all over the world (see 
e.g. http://www.worldwidewebsize.com/). This number 
starts to become even larger than 10 10 neurons in the 
brain. Each neuron can be viewed as an independent 
processing unit connected with about 10 4 other neurons 
by synaptic links (see e.g. [2-4]). About 20% of these 
links are unidirectional [5] and hence the brain can be 
viewed as a directed network of neuron links. At present, 
more and more experimental information about neurons 
and their links becomes available and the investigation 
of properties of neuronal networks attracts an active in- 
terest of many groups (see e.g. [6-13]. 

The WWW is also a directed network where a site 
j points to a site i but no necessary vice versa. The 
classification of web sites and information retrieval from 
such an enormous data base as the WWW becomes a 
formidable challenge of modern society where the search 
engines like Google are used by internet users in everyday 
life. An efficient way to classify and extract the informa- 
tion from WWW is based on the PageRank Algorithm 
(PRA), proposed by Brin and Page in 1998 [14], which 
forms the basis of the Google search engine. The PRA 
is based on the construction of the Google matrix which 
can be written as (see e.g. [15] for details): 

G = aS + (1 - a)~E/N . (1) 

Here the matrix S is constructed from the adjacency ma- 



World Wide Web; trix A of directed network links between N nodes so that 
Sij = Aij j Akj and the elements of columns with only 
zero elements are replaced by l/N. The second term in 
r.h.s. of (1) describes a finite probability 1 — a for WWW 
surfer to jump at random to any node so that the matrix 
elements Eij = 1. This term with the Google damping 
factor a stabilizes the convergence of PRA introducing 
a gap between the maximal eigenvalue A = 1 and other 
eigenvalues A^. As a result the first eigenvalue has Ai = 1 
and the second one has | A2 1 < a. Usually the Google 
search uses the value a = 0.85 [15]. By the construction 
J2i Gij — 1 so that the asymmetric matrix G belongs to 
the class of Perron- Frobenium operators which naturally 
appear in the ergodic theory [16] and dynamical systems 
with Hamiltonian or dissipative dynamics [17]. 

The right eigenvector at A = 1 is the PageRank vector 
with positive elements pj and ^2jPj — 1- The classifica- 
tion of nodes in the decreasing order of pj values is used 
to classify importance of WWW nodes as it is described 
in more detail in [15]. The PageRank can be efficiently 
obtained by a multiplication of a random vector by G 
which is of low cost since in average there are only about 
ten nonzero elements in a typical line of G of WWW. 
This procedure converges rapidly to the PageRank. 

Fundamental investigations of the PageRank proper- 
ties of the WWW have been performed in the com- 
puter science community (see e.g. [18-23]; involvement 
of physicists is visible, e.g. [24], but less pronounced). It 
was established that the PageRank is satisfactory charac- 
terized by an algebraic decay pj ~ 1/j^ with j being the 
ordering index and /3 « 0.9; the number of nodes with 
the PageRank p scales as N n ~ 1 jp v with the numerical 
value of the exponent v = 1 + 1/(5 rj 2.1 [15, 18]. It is 
known that such type of algebraic dependencies appear 
in various types of scale- free networks [25]. The PageR- 
ank classification finds its applications not only for the 
WWW but also for the network of article citations in 
Physical Review as it is described in [26, 27]. This shows 
that the approach based on the Google matrix can be 
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FIG. 1: Distribution of ingoing (left panels) and outgoing 
(right panels) links k: Pi„ and Pout give number of nodes 
with n ingoing and outgoing links respectively. Top panels: 
unweighted links; bottom panels: weighted links. 



applied to vary different types of networks. 

In this work we construct the Google matrix G for a 
model of brain analyzed in [11]. The properties of the 
spectrum and the eigenstates of G are described in the 
next Section II. The results are discussed in Section III. 



II NUMERICAL RESULTS 

To construct the Google matrix of brain we use a di- 
rected network of links between N = 10 4 neurons [28] 
generated from the brain model [11]. In total there are 
Ni = 1960108 links in the network. They form N out out- 
going links and AT in ingoing links (Ni — N out — Ni n ), so 
that there are about 200 outlinks (or ingoing) per neu- 
ron. These numbers include multiple links between cer- 
tain pairs of neurons; certain neurons have also links to 
themselves (there is one neuron linked only to itself). The 
number of weighted symmetric links is approximately 
9.8%. Due to existence of multiple links between the 
same neurons we constructed two G matrices based on 
unweighted and weighted counting of links. In the first 
case all links from neuron j to neuron i are counted as 
one link, in the second case the weight of the link is pro- 
portional to the number of links from j to i. In both 
cases the sum of elements in one column is normalized 
to unity. The distributions of ingoing (Pi n ) and outgoing 
(Pout) links are shown in Fig. 1. The weighted distribu- 
tion of ingoing links have a pronounced peaked structure 
corresponding to different regions of brain model consid- 
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FIG. 2: (Color online) PageRank pj for the Google matrix of 
brain model at a — 0.6, 0.85, 0.9, 0.95 and 0.99 shown by red, 
magenta green, blue and black solid curves (full curves from 
bottom to top at log 10 j — 0.3). The dotted black curve corre- 
sponds to a = 0.999 and demonstrates strong dependence of 
the PageRank on a in the vicinity of a = 1. Panels (a) and (b) 
correspond to unweighted and weighted links. For panels (a) 
and (b) the values of PAR are £ = 8223. and 8314., 6295. and 
6040., 5570. and 5046., 3283. and 3367., 28.4, 90.0, 1.09 and 
1.19 for a = 0.6, 0.85, 0.9, 0.95, 0.99, 0.999 respectively. Panels 
(c) and (d) show the dependence of the influence-PageRank 
p* (j) on j for the same values of a as for top panels respec- 
tively for unweighted and weighted links (for a > 0.6 there is 
a strong overlap of curves). 



ered in [11]. We note that the distribution of links is not 
of free-scale type. 

The dependence of the PageRank on a is shown in 
Fig. 2. For a = 0.999 almost all probability pj is concen- 
trated on one neuron. This is the only one neuron which 
is linked only to itself. With the increase of a up to 0.99 
the main part of probability is concentrated mainly on 
about 10 neurons that approximately corresponds to the 
number of peaks in the distribution of weighted ingoing 
links in Fig. 1 (bottom left panel). At the same time 
the PageRank has a long tail at large j where the prob- 
ability pj is practically homogeneous. For a = 0.6 the 
peak of probability at 1 < j < 10 is washed out and the 
PageRank becomes completely delocalized. We note that 
a derealization of the PageRank with a appears in the 
Ulam networks describing dynamical systems with dissi- 
pation [29, 30]. At the same time the WWW networks 
remain stable in respect to variation of a as it is discussed 
in [23, 31]. 

Recently, for the studies of procedure call network of 
the Linux Kernel [32] it was proposed to study the prop- 
erties of the importance-PageRank p*(J) which is given 
by the eigenvector at A = 1 for the Google matrix con- 
structed from the inverted links of the original adjacency 
matrix. It was argued that p*(j) can give an additional 



0.2 
0.1 



£ 0.0- 
-0.1- 



-0.2 




0.2- 
0. 1 - 



£ 0.0- • 
-0.1- 



-0.25 0.0 0.25 0.5 0.75 1.0 

ReX 



-0.2- 



-0.25 0.0 0.25 0.5 0.75 1.0 

ReX 




0.1>r 



8 0.0 




-0.1 4^ 

0.4 0.6 0.8 1.0 

ReX 

FIG. 3: (Color online) Spectrum of eigenvalues of the Google 
matrix G of brain at a — 0.99 in the complex plain A for (a) 
unweighted and (b) weighted links in the neuronal network. 
Panels (c) and (d) show zooms of data of panel (c). The color 
shows the degree of localization of eigenvectors of G being 
proportional to the value of PAR £ and changing from one 
(red/light gray) to maximal value (dark green/black). 



information about certain important nodes. Our results 
for p*(j) are shown in panels (c,d) of Fig. 2. They show 
that p*(j) is practically delocalized and fiat for all used 
values of a. This indicates that all nodes have practically 
equal importance. The popularity-importance correlator 
introduced in [32] and defined as k — N ^2 i p{i)p* {i) — 1 
is rather small (n m —0.009, —0.017 at a = 0.6 and 
k ~ —0.054, —0.065 at a = 0.85 for unweighted, weighted 
links respectively). This shows that there are no correla- 
tions between p and p* in our neuronal network that is 
similar to the Linux Kernel case. 

The spectrum Ai and the right eigenvectors ipi of the 
Google matrix of brain are defined by the equation 
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The spectrum of A is complex and is shown in Fig. 3. 
The color of points is chosen to be proportional 
to the PArticipation Ratio (PAR) defined as £ = 
(Erj I^O')! 2 ) 2 / Yjj \A{j)\ 4 - This quantity determines an 
effective number of sites populated by an eigenstate ipi, 
it is often used to characterize localization-delocalization 
transition in quantum solid-state systems with disorder 
(see e.g. [33]). The spectrum has eigenvalues with |Aj| 
being close to unity so that there is no gap in the spec- 
trum of A in the vicinity of A — 1 (we remind that the 
second term in the r.h.s. of (1) transfers Ai to cvAi keep- 
ing only one Ai = 1 [15]). This is different from the 
spectrum of random scale-free networks characterized by 
a large gap in the spectrum of A [34] . 

Compared to the spectra of the university WWW net- 
works studied in [31] the spectrum of G in Fig. 3 is more 
flat being significantly compressed to the real axis. In 




FIG. 4: Dependence of the density of states dW/d'y of G on 
the relaxation rate 7 for unweighted (pluses) and weighted 
(circles) links in the neuronal network. 





FIG. 5: Dependence of PAR £ on relaxation rate 7 at a = 0.85 
for (a) unweighted and (b) weighted links in the neuronal 
network. 



this respect our neuronal network has certain similarity 
with the spectra of vocabulary networks analyzed in [31] 
(see Fig. 1 there). At the same time the spectrum of 
G matrix of brain has visible structures in the eigen- 
values distribution in the complex plane of A while the 
vocabulary networks are characterized by structureless 
spectrum. The spectrum of Fig. 3 has global properties 
being similar to those of the Ulam networks considered 
in [29]. It is interesting to note that the spectra of un- 
weighted and weighted networks of brain have similar 
structure. This supports the view of structural stability 
of the spectrum of G matrix. 

It is useful to determine the relaxation rate of eiges- 
tates by the relation 7 = — 21n|A|. The dependence of 
density of states dW/dry on 7 is shown in Fig. 4 (the den- 
sity is normalized to unity so that J dW/ d'yd'j = 1 cor- 
responds to N — 10 4 states). The distribution in 7 has a 
pronounced peak at 7 rs 5, the density of states at small 
7 < 1 is relatively small (this is also seeing in Fig. 3). 
The comparison of unweighted and weighted links shows 
the stability of the density distribution in respect to such 
modification of links. 

The dependence of the PAR £ on 7 is shown in Fig. 5 
(we note that except of the PageRank £ is independent 
of a due to the unity rank of matrix E, see e.g. [15, 29]). 
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FIG. 6: Dependence of PAR £ of the PageRank on parameter 
a for (a) unweighted and (b) weighted links in the neuronal 
network. 
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FIG. 7: Distribution of PageRank values p and p* for all sites 
for (a) unweighted and (b) weighted links in the neuronal 
network at a — 0.85. 



The PageRank value of £ at 7 = is very large being more 
than half of the total number of neurons N = 10 4 . It is 
clear that this corresponds to a delocalized state. The 
eigenstates with < 7 < 2 have relatively small £ < 10 3 
being close to a localized domain while eigenstates with 
2 < 7 < 10 have £ > 10 3 being delocalized on the main 
part of the network; the states with 7 > 10 enter in 
the localized domain. For a > 0.99 the PAR is close 
to £ ps 1. Taking as a criterion that the derealization 
takes place when £ > N/2 we obtain that the PageRank 
becomes delocalized at a c « 0.9 (see data of Figs. 2,6). 
The global dependence of the PAR £ of the PageRank on 
parameter a is shown in Fig. 6 with a sharp derealization 
of £ for a < a c . Of course, the above analysis should be 
considered as an approximate one since the localization 
properties should be studied in dependence on the system 
size N while we consider only one size of N. 

Finally, following the approach proposed in [32], we 
show in Fig. 7 the distribution of PageRank values p and 
p* for all sites. Such kind of distributions can be rather 
useful in determining sites which have maximal values of 
p and p* at the same time. However, a detailed analysis of 
the properties of this distribution would require networks 
with a larger size N where statistical fluctuations are 
smaller. 



Ill DISCUSSION 

In this work we studied the properties of the Google 
matrix of a neuronal network of the brain model dis- 
cussed in [11]. For this network of 10 4 neurons we found 
that the spectrum of the Google matrix has a gapless 
spectrum at a = 1 demonstrating certain similarities 
with the spectra of university WWW networks and vo- 
cabulary networks studied in [31]. At the same time our 
neuronal network shows signs of derealization transition 
of the PageRank at the Google damping factor a c w 0.9 
which was absent in the networks studied in [23, 31]. A 
similar transition in a was detected in the Ulam net- 
works generated by dissipative dynamical maps [29] . We 
attribute the appearance of such derealization transition 
to a large number of links per neuron (200) which is by 
factor 10 larger than in the WWW networks (20). 

Of course, our studies have certain limitations since 
we considered only a fixed size neuronal network and 
since this network is taken from a model system of brain 
analysed in [11]. Another weak point is that we do not 
consider the dynamical properties of the network which 
are probably more important for practical applications. 
Nevertheless, the spectral properties of G matrix can be 
rather useful. Indeed, the gapless spectrum of A shows 
that long living excitations can exist in our neuronal net- 
work. Such relaxation modes with small rates 7 can be 
the origin of long living oscillations found in numerical 
simulations [11]. It is quite possible that the properties 
of spectra of G can help to understand in a better way 
rapid relaxation processes and those with long relaxation 
times. We conjecture that the rapid relaxation modes 
correspond to relaxation of local groups of neurons while 
long living modes can represent relaxation of collective 
modes representing dynamics of human thoughts. The 
dynamics of such collective modes can contain significant 
elements of chaotic dynamics as it was discussed in the 
frame of the concept of creating chaos in [35] . 

It is possible that the brain effectively implements dy- 
namics described by the evolution equation dip/dt — Gip 
which without perturbations converges to the steady- 
state described by the PageRank (which may be linked 
with a sleeping phase). External perturbations give ex- 
citations of other eigenmodes of G discussed here. The 
evolution of these excitations will be significantly affected 
by the spectrum of G. 

Further development of the Google matrix approach to 
the brain looks to us to be rather promising. For exam- 
ple, a detection of isolated communities and personalized 
PageRank, represented by other types of matrix E in (1), 
is under active investigation in the computer science com- 
munity (see e.g. [15, 23]). Such type of problems can find 
their applications for detection of specific quasi-isolated 
neuronal networks of brain. The usage of real neuronal 
networks, similar to those studied in [6-10, 13], in com- 
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bination with the Google matrix approach can allow to 
discover new properties of processes in the brain. The de- 
velopment of parallels between the WWW and neuronal 
networks will give new progress of the ideas of John von 
Neumann. 

We thank E.M. Izhikevich for providing us with the 
data set of links between neurons [28] in the brain model 
[11]. 
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